Disorder Effect on Spin Excitation in Double Exchange Systems 
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Spin excitation spectrum is studied in the double exchange model in the presence of disorder. Spin 
wave approximation is applied in the lowest order of 1/S expansion. The disorder causes anomalies 
in the spin excitation spectrum such as broadening, branching, anticrossing with gap opening. The 
origin of the anomalies is the Friedel oscillation, in which the perfectly polarized electrons form the 
charge density wave to screen the disorder effect. Near the zone center q = 0, the linewidth has 
a q linear component while the excitation energy scales to q 2 , which indicates that the magnon 
excitation is incoherent. As q increases, there appears a crossover from this incoherent behavior to 
the marginally coherent one in which both the linewidth and the excitation energy are proportional 
to q 2 . The results are compared with experimental results in colossal magnetoresistance manganese 
oxides. Quantitative comparison of the linewidth suggests that spatially-correlated or mesoscopic- 
scale disorder is more relevant in real compounds than local or atomic-scale disorder. Comparison 
with other theoretical proposals is also discussed. Experimental tests are proposed for the relevance 
of disorder. 

PACS numbers: PACS numbers: 75.30.-m, 75.47.Lx, 75.47.Gk 



I. INTRODUCTION 

Since rediscovery of the colossal magnetoresistance 
(CMR) phenomena, perovskite manganese oxides have 
attracted much attentioni2iiMi£i£i£ This family of com- 
pounds shows a variety of physical properties in magnetic 
(ferromagnetism or antiferromagnetism) , charge (metal, 
insulator, or charge ordering), and orbital/lattice (orbital 
ordering and/or Jahn- Teller distortion) degrees of free- 
dom. Phase transitions between different ordered phases 
can be controlled by chemical substitution of the A site 
in the general chemical formula j4Mn03. The magnitude 
of the CMR effect also changes drastically for different A- 
site ions, hence, one of the most important issues in this 
field is to understand the role of the ^4-site substitution. 

The ferromagnetic metallic state, which is widely sta- 
ble at low temperatures in these compounds,^ is basi- 
cally understood by the Zener's double exchange (DE) 
mechanismii2iiiii£ The kinetics of itinerant electrons in 
the Mn e g bands strongly correlates with localized spins 
in the t2 g levels through the large Hund's-rule coupling; 
ferromagnetism of the localized spins leads to high con- 
ductivity of the electrons, and vice versa. Thus, the issue 
is what is an additional element to the DE interaction 
which is necessary to explain deviations from this canon- 
ical DE behavior and instabilities towards a variety of 
different phases in the ^4-site substituted materials. 

Among various manganese oxides, it has been recog- 
nized that Lai-^Sr^MnOs (LSMO) at x ~ 0.3 is a canon- 
ical DE system. 13 Namely, the simple DE model can ex- 
plain magnetic, transport and optical properties of this 
material quantitatively. For instance, the experimental 



value of the Curie temperature Tq is well reproduced 
by recent theoretical calculations in the DE model^^ 
This material shows relatively high Tc compared to other 
manganese oxides, which also suggests the stability of 
ferromagnetic metal due to the primary role of the DE 
mechanism and the irrelevance of other additional fac- 
tors. This canonical DE system gives a good starting 
point to examine effects of the ^4-site substitution. 

There are two different ways of the systematic ^4-site 
substitution. One is the substitution with different ionic 
valences, such as the control of x in Lai-xSr^MnOa. The 
mixture of trivalent (La 3+ ) and divalent (Sr 2+ ) ions effec- 
tively changes the nominal valence of manganese ion as 
M n (3+i)+ i which controls the electron density in the Mn- 
O network. The other substitution is by different ionic 
radii at a fixed valence, such as Ai_ x A'' x MjiOz for dif- 
ferent combinations of, for instance, ^4=La,Pr,Nd,Y and 
j4'=Ba,Sr,Ca at a fixed x. These two different substi- 
tutions cause systematic and drastic deviations from the 
canonical DE behavior in LSMO compounds at x ~ 0.3. 

We discuss mainly the latter ^4-site substitution with 
the ionic-radius control in the following. The change of 
the ionic radii is a sort of the chemical pressure, which 
modifies length and angle of Mn-O-Mn bonds. This leads 
to the change of effective transfer integrals between Mn 
ions, namely, the bandwidth of the electrons. Thus, this 
substitution has been often called the bandwidth control. 
This terminology 'bandwidth control', however, might be 
misleading: There are many experimental facts which 
cannot be explained only by the change of the bandwidth. 
One is a rapid decrease of Tq compared to the bandwidth 
changei 16 : 17 Tq decreases as the averaged radius of the A 
ions decreases; however, the change of Tq is much larger 
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than that of the bandwidth which is estimated from the 
structural change. For instance, from Lao.7Sr .3Mn03 to 
La . 7 Cao.3Mn03 (LCMO), T c decreases by about 30% 
while the estimated bandwidth decreases by less than 
2%iA£ Since, in the DE theory, Tc is proportional to the 
bandwidth, the large change of Tc cannot be explained 
by the bandwidth change alone. There should be a hid- 
den parameter which strongly suppresses the kinetics of 
the electrons. 

In this work, we focus on the quenched disorder for 
a candidate of the hidden parameter in this ^4-site sub- 
stitution. In general, oxides are known to be far from 
perfect crystals. Moreover, in these manganese oxides, 
the disorder is inevitably introduced since they are solid 
solutions of different j4-site ions. 

The importance of the disorder has been pointed out 
in several experimental results: (i) Tc changes not only 
as the average of the ionic radius but also as its standard 
deviation^ Even if the average is the same, Tc becomes 
lower for larger standard deviation, (ii) The residual 
resistivity, namely, the value of the resistivity extrapo- 
lated to zero temperature, becomes larger for lower- Tc 
compoundsii2i2S (iii) ^4-site ordered materials, in which 
different A ions form a periodic ordered structure, exhibit 
higher Tc compared to compounds with random distri- 
bution of A ions in the same chemical formulai^iiSi^ 
All these experimental results suggest that the disorder 
in the random mixture of different size ions scatters the 
itinerant electrons and suppresses the kinetics of them. 

In this paper, we discuss the disorder effect on spin dy- 
namics in the ^4-site substituted manganites. The spin 
dynamics is another important indication of a hidden 
parameter. Spin excitation spectrum shows qualitative 
changes for the ^4-site substitution. In compounds with 
relatively high Tc such as LSMO, the spin excitation 
shows a cosine type dispersion, which is similar to that 
of Heisenberg spin systemsi^2^ This behavior is well de- 
scribed by the DE mechanism alone^ On the contrary, 
in compounds with low Tc such as LCMO, the spec- 
trum shows significant deviations from this form, e.g., 
some anomalies such as broadening, softening, and gap 
openingi^£i2£i22iSi To explain these anomalies is one of 
the crucial test for elucidating the hidden parameter. 

We will demonstrate here that the disorder gives a 
comprehensive understanding of systematic changes of 
the spin excitation spectrum as well as other experimen- 
tal results described above. Several other mechanisms 
have been proposed, for instance, antiferromagnetic in- 
teractions between localized spins^ orbital degrees of 
freedom in the e g bands^the electron-lattice coupling, 31 
and the electron-electron correlationi^i 3 ^^ Through the 
detailed analysis of the spectrum, we will show that the 
disorder appears to be promising among these scenarios. 
A part of the results has been published in the previ- 
ous publicationsi^LSiiS More systematic and extensive 
analyses are presented in this paper. 

This paper is organized as follows. In Sec. [HJ we de- 
scribe methods to study spin dynamics in the DE model. 



To be self-contained, we derive and summarize some an- 
alytical expressions which have been developed in the 
previous publicationsi^LSiiS Numerical results are pre- 
sented in Sec. IIIII for wide regions of parameters. We 
find remarkable anomalies in the spin excitation spec- 
tra and discuss their origin in detail. Comparison with 
the analytical results is also examined. The results are 
compared with experimental results as well as with other 
theoretical ones in Sec. IIVI We propose possible experi- 
ments to test the relevance of the disorder. Section M is 
devoted to summary and concluding remarks. 



II. FORMULATION 
A. Model 

In this work, we study effects of the disorder on the 
spin dynamics of the DE modelii2*ii Our Hamiltonian is 
given in the form 

i<Cj (7— ± i 
ia 

Here, the first term describes the electron hopping be- 
tween i and jth site, the second term represents the 
Hund's-rule coupling between the electron spin a (vector 
of Pauli matrix) and the localized spin S, and the last 
term contains the on-site potential Si. 

In model Q , we take account of the quenched disorder 
in two ways. One is the off-diagonal disorder in the trans- 
fer integral tjj and the other is the diagonal disorder in 
the potential energy e^. Hereafter, we call the former the 
bond disorder and the latter the on-site disorder, respec- 
tively. The following discussions in this Sec. [H] do not 
depend on the detailed functional form of the distribu- 
tion of the random variables fcy and E{ . The distribution 
will be specified in Sec. IIII Bl for numerical calculations 
in Sec.lTTTl 



B. 1/5 expansion 

We apply the spin wave approximation in the lowest 
order of 1/5 expansion to model (S is the mag- 
nitude of the localized spin S). We consider a fer- 
romagnetic ground state with perfect spin polarization 
S t = (0, 0, S)^ and calculate the one-magnon excitation 
spectrum. In the absence of disorder, the formulation is 
given in Ref. |2(| Here, we extend the method to the case 
with disorder. 

In the presence of disorder, it is convenient to 
work with the real-space representation instead of the 
momentum-space one. First, we explicitly diagonalize 
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FIG. 1: Spin wave self-energy in the lowest order of the 1/5 
expansion. See Ref. [2(J for details. 



the Hamiltonian matrix of model JQ) for a given config- 
uration of the disorder {tij,Ei}. We have 

E Hi] ({tjj , £»})¥W (j) = {En - <TJ-R)ipn<j(i), (2) 



where ip ntT is the nth eigenstate with the eigenenergy E n . 
Electron Green's function is given by 



uj - (E n - aJft - fx) + irj sgnw ' 



(3) 



where fi is the chemical potential and 77 is an infinitesi- 
mal for convergence. In the lowest order of the 1/5 ex- 
pansion, following Ref. [2^, the magnon self-energy is ob- 
tained from the electron spin polarization function shown 
in Fig. ^ The contribution from Fig. ^ (a) is given by 



s 



(4) 



where /„+ is the fermi distribution function for up-spin 
states. The other contribution from Fig. ^ (b) is given 
by 

nxy, \ 2 ^ (p n+ {i)<p* n+ (j)(p m+ (j)ip; n+ (i) 



w + E n — E m — 2 Jh 



(5) 

In the limit of Jh 3> t, which is realistic in CMR man- 
ganites, the summation of Eqs. J3J and (JSJ ends up with 

mn 

x(E m - E n - uj). (6) 
Here, we assume the orthonormal relation 

m 

and drop the spin index for simplicity. 

Magnon Green's function is calculated by this self- 
energy in the form 



Dij(u>) = [ oj - IIjj-(u>) +ir) ] 1 



(8) 



Spin wave excitations are obtained from the poles of the 
Green's function JSJ. Since the self-energy in Eq. © 



is proportional to 1/5, the positions of the poles are at 
uj ~ 0(1/ S). Hence, within the lowest order of 1/5 ex- 
pansion, spin wave excitation is determined by the static 
part of the self-energy Hy (uj — 0) as discussed in Ref.Erj 
The matrix inversion in Eq. (|SJ) is easily obtained by using 
the eigenvectors ipi and the eigenvalues uji which satisfy 



^Hy(w = 0) • ipi(j) = ujitpi(i). 



(9) 



Finally, the spectral function for a given configuration 
of {tij is calculated as 

A <*> u ) = -^E~ im ^>) el£Hri ~ rj3 
= ^ElE^oy^lV-wo, (10) 

where N is the system size. The spin excitation spectrum 
is obtained by averaging A(q, uj) for random configura- 
tions of {tij, Si}. In the same manner, all the physical 
quantities are to be averaged for random configurations 
although we do not explicitly remark in the following dis- 
cussions. 

To summarize, the spin excitation spectrum in the DE 
model is calculated up to the lowest order of the 1/5 ex- 
pansion as follows, (i) Diagonalize the Hamiltonian 
for a configuration of the disorder {tij,£i}, and obtain 
the electron wavefunction (p n and eigenenergy E n . (ii) 
Calculate the magnon self-energy © by using <p n and 
E n . (iii) Diagonalize the self-energy as in Eq. JJjJ, and 
obtain the eigenvalues u>i and eigenvectors ipi. (iv) Cal- 
culate A(q,uj) according to Eq. (|ll)fl . (v) Repeat (i)-(iv) 
for different random configurations and take the random 
average of A(q, uj). 

C. Analytical formula 

In this section, we derive some analytical expressions 
for the spin wave excitation. The following framework is 
general and does not depend on the details of the disorder 
and the spatial dimension of the system. The derived 
formulas are useful to discuss the numerical data in the 
following sections. 

1. Magnon self-energy 

First, we examine the static part of the magnon self- 
energy. We can rewrite Hy(u; = 0) in the form 



(°) - go ~ % E /> ' / < 



(ii) 



where Bij = J2 n /nVn(i)VnW- H ere : we use the follow- 
ing relations 

E m ip m (i)(f* m (j) = ^2n ik <Pm(k)ip* m (j) 

m mk 
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= Hij, (12) 

n nk 

= J2 H 3kB ki , (13) 

k 

with the eigenequation (J2Jl and the orthonormality Q). 

From Eq. ^1 If) , we can easily show that the self-energy 
satisfies the sum rule in the form 

E n « (°) = ^ (E n -> ,! - - E = °- ( 14 ) 

3 3 k 

Thus, the eigenvalue in Eq. 10 always becomes zero at 
q = even in the presence of disorder. This indicates 
that there is a gapless excitation at q = 0. This is con- 
sistent with the Goldstone theorem. 

We note also that the matrix element of Ily (0) is given 
by the transfer energy of the electrons. By using the 
relation 

tUj B 3i = E fn'Pli^ijfnU) = ( 15 ) 
n 

the matrix elements are written as 

II i7 y(0) = -ty( c Jc,-)/25 = -2SJij, (16) 
n«(0) = 25 ^ J i3 -, (17) 

where we define 

J ij =t i Mci)/*S a - ( 18 ) 

The bracket represents the expectation value in the 
ground state for a given configuration of the disorder. 
We will show in Sec. Ill Dl that is the exchange cou- 
pling of the corresponding Heisenberg spin model. Equa- 
tions (|16fl and (|17fl indicate that the following summa- 
tions give the kinetic energy of electrons; 

E n «(°) = -E n »(°) = <^kin)/25 (19) 

where 7ikin is the first term in the Hamiltonian Q). 

From Eqs. (|16|l and Q17[l. we note that the magnon 
self-energy at cj = is real in general since (cjcj) is real 
for the fully-polarized ferromagnetic ground state with- 
out degeneracy. Hence, as seen in Eqs. © and i fTUj l. up 
to 0(1/5), A(q,u>) for a given configuration of disorder 
describes a well-defined quasi-particle excitation with in- 
finite lifetime (zero linewidth) . Lifetime due to the decay 
of the quasi particle is obtained from the imaginary part 
of the poles of Eq. ©, however, it is a higher order term 
in the 1/5 expansion. (This higher order correction will 
be discussed in Sec. IIVB 21 ) On the other hand, in the 
presence of disorder, we have to take the random average 
for different realizations of disorder configurations. Since 
different configurations give different A(q,Lu), we obtain 



a distribution of the excitation spectra after the random 
average. The distribution gives rise to a finite linewidth, 
which can be estimated, for instance, by the standard de- 
viation of the excitation energy. This linewidth is due to 
the so-called inhomogeneous broadening. In the presence 
of the disorder, this broadening is the lowest order effect 
in the 1/5 expansion. We will analyze this linewidth in 
the following sections. 

2. Spectral function analysis 

If the spin excitation spectrum is single-peaked and 
does not show any splitting, we can apply the follow- 
ing spectral function analysis. We will show in Sec. IIIII 
that this is the case near the zone center q — 0. Let us 
consider the mth moment of the spectral function as 

/>oo 

n< m) = / ui m A(q,u) duj. (20) 
Jo 

Using Eqs. 10 and (|TU)l . we obtain 

n «° = jf E E • • • r r • 

ij k\k2---k rn — \ 

(21) 

Thus the moment Qg 1, is a Fourier transform of the 
magnon self-energy (IIy(0)) m . 

The first moment in Eq. (|20|l gives the averaged exci- 
tation energy as 

w*w{q) = n§ ) . (22) 

As mentioned in the last part of Sec. Ill Cll the linewidth 
is estimated by the standard deviation of the excitation 
energy which is given by the second moment in the form 

7(«) = [n«-(nJJ) 9 ] 1/a . (23) 

3. Excitation energy 

Here, we analyze the excitation energy of Eq. Q22p. 
Substituting Eqs. (2J and JUJ into Eq. we obtain 

- S w( Q ) = ^E n ^ ( " 3) - ( 24 ) 

When model JTJ has the electron hopping only between 
the nearest neighbor sites, the magnon self-energy 
has nonzero matrix elements only for i — j and i = j + rj 
as shown by Eqs. (|16fl and i|17|) . Here r\ is a displacement 
vector to the nearest neighbor site. Then the excitation 
energy is given as 

^sw(g) = iE( n » + E n ^ eiq ' V ) 

v 

= E n(77) (e'^-l), (25) 
V 
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where we define the site-averaged quantity II(tj) = 
Tli,i+rj /N and use the sum rule Eq. (fT^l) . Even in 
the presence of the disorder, the symmetry for the direc- 
tion of 77 is expected to be recovered after the random 
average, therefore we assume 



n(i,) = i£ fl fa) = - A . 
z v 



(26) 



where z is the number of coordinates. Finally, we obtain 
the following form 



Wsw (q) = A^(l-e i 9- 7 ?). 
V 



(27) 



This indicates that on hypercubic lattices the excitation 
spectrum shows the cosine dispersion even in the presence 
of disorder. In the limit of q = \q\ — > 0, we have 



UJ sw {q) ~ Aq 2 



(28) 



Thus, the quantity A gives the spin wave stiffness. In 
our analysis, the stiffness is proportional to the kinetic 
energy of electrons as 



A : 



1 



1 



(H 



kin/ 



v 



2SzN'- 



(29) 



where we use Eq. lfTU|) . 



4- Linewidth 

Next, we analyze the linewidth of Eq. I|23|l . Similarly 
to the derivation of Eq. I)25H , the second moment of the 
spectral function is written in the form 

^ 2) =E Mm^il-^^il-e^), (30) 



where 



n 2 (i7i,Tj 2 ) = ^X!n i+T?1 ,in M+ r7 2 , (3i) 



In the limit of q — > 0, we have 

7 2 (q) ~cig 2 + c 2 g 4 , 

where 

d = 4 (5A 2 , 

4 , 
c 2 = A 2 - - <5A 2 - A^ 



(34) 

(35) 
(36) 



Therefore, in the presence of disorder, the linewidth 
shows a q linear behavior in the small q regime. 



5. Incoherence of spin wave excitation 

From the above analysis, we obtain lu sw cx q 2 and 7 oc q 
in the limit of q — > as in Eqs. (J2EJ and l|3"4"j). Thus, the 
linewidth 7 becomes larger than the excitation energy 
w sw at q ~ 0, which indicates that the spin wave excita- 
tion becomes incoherent or localized. 

This incoherent behavior comes from local fluctua- 
tions. The q linear coefficient in 7 in Eq. Q35J1 is given by 
the difference of Il 2 for different directions of r/, which is 
proportional to 

ci oc n 2 (r7,rj) - n 2 (rj, -77) 

oc y^(u i+11t t - ii^-Tj) 2 



(37) 



where we use Eq. I|16|l . From the definition of Jy in 
Eq. iJTSJl, Eq. I|37|l shows that the q linear coefficient C\ 
is proportional to the local fluctuations of the transfer 
energy of electrons. 

As q becomes large, we have a crossover from this in- 

1/2 

coherent regime where 7 cx w S w oc q to the marginally 
coherent regime where 7 cx ui sw cx q 2 . This crossover is 
determined by the coefficients c\ and c 2 in Eqs. I|35|l and 
(|36|l . These coefficients depend on the details of parame- 
ters in the system, especially on the type of the disorder, 
as will be discussed in Sec lIIIE2l 



and r) x and T] 2 are the displacement vectors to the near- 
est neighbor site. In the absence of disorder, we have 
i+r].i = ^i,i+r] = — A. Hence, the linewidth 7 becomes 



n 



>( 2 ) 



q _ (n { q } ) 2 . In the presence of disorder, the 



zero since f2 

linewidth becomes finite. For simplicity, we consider here 
the linewidth in a special direction such as q = (q, 0, 0, • • •) 
on a hypercubic lattice. Since fl 2 in Eq. (|31l) generally 
depends on the relative direction between rj 1 and ry 2 , we 
represent two different cases as 



Tl 2 {r] 1 ,r] 2 ) = A 2 ±SA 2 for r) l =±r} 2 . 
Then, Eq. I|3U|) is rewritten as 



n 



(2) 



4( A 2 (l - cosq) 2 + 5A 2 sin 2 q ). 



(32) 



(33) 



D. Correspondence to Heisenberg model 

Here, we discuss the relation between the spin exci- 
tation spectra in the DE model and those in the 
Heisenberg spin model whose Hamiltonian is given by 
^Hcis = — 2X)i<j JijSi ■ Sj. As explicitly shown in Ap- 
pendix ^ the static part of the magnon self-energy 
corresponds to the effective spin-wave Hamiltonian for 
the Heisenberg model within the leading order of the 1/S 
expansion. There, Eq. (|18l) gives the relation between the 
kinetics of itinerant electrons in the DE model (JIJ and 
the exchange coupling in the Heisenberg model. This 
relation provides us physical intuitions to understand our 
results. 
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In the absence of disorder, when we consider the trans- 
fer fcy = t only for the nearest neighbor sites, Jy in 
Eq. (|18|l becomes a constant J which is nonzero only for 
i = j + V- Thus, as noted in Ref. l2(i the spin excitation 
spectrum in the pure DE model with infinite Jh is equiv- 
alent to that in the Heisenberg model with the nearest 
neighbor exchange Jy = J, and is given by the simple 
cosine dispersion. 

Disorder introduces some randomness in the expecta- 
tion value of (c\cj) as well as the transfer integral iy. 
This leads to the nonuniform exchange coupling Jy in 
the corresponding Heisenberg model. It is instructive to 
compare the spectrum for the DE model and that for 
the Heisenberg model with random exchange couplings 
Jij which no longer satisfy the relation l|18l) . (See also 
Sec. IIII D3l l For the random Heisenberg model, the spin 
wave stiffness A is given by the random average of the 
exchange constant Jy. Moreover, the q linear term in 
the linewidth is determined by local fluctuations of the 
random exchange Jy. Thus, these aspects are commonly 
observed both in the random DE model and in the ran- 
dom Heisenberg model even if the relation l|18fl does not 
hold between them. In other words, in the long wave- 
length limit of q — > 0, the magnetic excitation has a sim- 
ilar structure between these two models, although one is 
the itinerant electron system and the other is the local- 
ized spin system. However, for moderate or large g, the 
excitation spectrum of the DE model in the presence of 
disorder shows remarkable differences from that of the 
random Heisenberg model as shown in Sec. II I II There, 
the fermionic properties of itinerant electrons play a key 
role through the expectation value (c\cj). 

III. RESULTS 

A. Overview 

In this section, we show numerical results which are 
obtained based on the method described in Sec. Ill Bl Af- 
ter the specification of the functional form of the disorder 
distribution in Sec. IIII Bl we show that the spin excita- 
tion spectrum exhibits some anomalies due to the dis- 
order in Sec. llHCl We examine the anomalies in detail 
by changing parameters such as strength and type of the 
disorder, spatial dimensions, and doping concentration. 
We discuss the origin of these anomalies in Sec. IIII Dl In 
Sec. IIII El we examine quantitative aspects of the excita- 
tion spectrum by applying the spectral function analysis 
in Sec. Ill C"2l to the numerical results. We compare them 
with the analytical expressions obtained in Sec. Ill CI 

B. Distribution form of disorder 

In the following sections, we consider the on-site dis- 
order Si = Ss and the bond disorder £y — t + St only for 
the nearest neighbor sites where 8s and St obey one of 



FIG. 2: Spin excitation spectra at x = 0.3 in the presence of 
the on-site disorder in the binary distribution with (a) A = 0, 
(b) A = 0.1, (c) A = 0.2, and (d) A = 0.3, respectively. 



the following three distribution functions. The first one 
is the binary distribution in the form 

P(x) = ±[5(A) + S(-A)], (38) 

namely, x = +A or —A in the equal probability. The 
second one is the Gaussian distribution in the form 




The last one is the box distribution in the form 

^) = ^(i+p)(i-e(i- P )), (40) 

where p — \/3A and Q(x) is the Heaviside step function. 
The normalizations are taken to give the same second 
moment as J x 2 P(x)dx — A 2 . 

In the following calculations, we consider the hyper- 
cubic lattice in d dimensions, and we take the half- 
bandwidth W = zt = 2dt for Jh = A = as an energy 
unit. We change the value of A as a parameter typi- 
cally up to 0.5VF. It is difficult to determine the realistic 
value of A in low Tq manganites. A rough estimate has 
been given by the first principle calculation, which shows 
that A becomes the same order of magnitude of the half- 
bandwidth WM> 



C. Anomaly in spin excitation spectrum 

1. Dependence on strength and type of disorder 

We first show the change of the spin excitation spec- 
trum by controlling the strength of the disorder in Fig. [5] 
Here, we fix the doping concentration at x = 0.3, where 
the hole density x is defined as x = 1 — (J2i c l°i/N) 
(we drop the spin index because the ground state is per- 
fectly polarized). Numerical calculations are performed 
for N = 20 x 20 x 20 site clusters under the periodic 
boundary conditions. The results are for the on-site dis- 
order whose distribution is given by the binary form of 
Eq. (|38|l . The spectrum is obtained by the random av- 
erage of A(q,u>) in Eq. (|10|l over 16 different realizations 
of random configurations for each value of A. The gray- 
scale contrast shows the intensity of the spectrum. 

In the case of A = 0, i.e., in the absence of the dis- 
order, the spin excitation spectrum is given by a cosine 
dispersion as shown in Fig. (a) . When we switch on the 
disorder, the excitation shows a finite linewidth which be- 
comes large as A increases. Moreover, as clearly seen in 
Figs. 121(c) and (d), the spectrum shows some anomalies 
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FIG. 3: Spin excitation spectra at x = 0.3 in the presence of 
the on-site disorder with A = 0.3 in (a) the Gaussian distri- 
bution and (b) the box distribution. 



FIG. 4: Spin excitation spectrum at x — 0.3 in the presence 
of the bond disorder in the binary distribution with A = 0.15. 

in the large q region. There, we have additional broad- 
ening as well as branching, that is, an emergence of an 
additional branch which has significantly lower energy 
than the original cosine like excitation. For large val- 
ues of A, the lower branch has a substantial weight as 
shown in Fig. 0(d). If one follows only the lower branch, 
the spectrum appears to show softening near the zone 
boundary. 

These anomalous features are commonly observed ir- 
respective of the type of the disorder. Figure shows the 
results for different distribution functions; (a) is for the 
Gaussian distribution of Eq. (|39|1 and (b) is for the box 
distribution of Eq. I|4U|). In both cases, we have similar 
anomalous features to those in Fig. [2J 

We also examine the case of the bond disorder in Fig. 21 
Here, we consider the binary distribution (|38[> for St. In 
this case also, the spectrum shows additional broadening 
which is more prominent and makes a branching obscure 
compared to the cases of the on-site disorder in Figs. [21 
and|31 This is related to strong incoherence in the case of 
the bond disorder which will be discussed quantitatively 
in Sec. iDIEl 



2. Spatial dimension dependence 

Next, we examine systems in less spatial dimensions, 
namely, in two and one dimensions. In manganites, 
there are compounds which have strong two-dimensional 
anisotropy, such as the bilayer materials ^M^O^ and 
the single-layer materials ^MnO^. Later, we will discuss 
the quantitative comparison between our results and ex- 
perimental ones in these compounds in Sec. IIVI Here, 
we calculate 40 x 40 site clusters for 50 different real- 
izations of random configurations in the two-dimensional 
(2D) case, and 256 site clusters for 500 configurations in 
the one-dimensional (ID) case. 

Figure shows the excitation spectra in 2D and ID 
systems at x = 0.3 for the on-site disorder in the bi- 
nary distribution of Eq. I|38|) with A = 0.3. In both 
dimensions, the spectra show the anomalies clearly as in 
three-dimensional (3D) systems in Sec. IIII C II The im- 
portant aspect is that in lower dimensions the anomalous 



FIG. 5: Spin excitation spectra in (a) two dimensions and (b) 
one dimension. The results are for x = 0.3 with the on-site 
disorder in the binary distribution with A = 0.3. 



FIG. 6: Spin excitation spectra in three dimensions for (a) 
x = 0.5, (b) x — 0.4, (c) x — 0.3, and (d) x — 0.2, respec- 
tively. The results are for the on-site disorder in the binary 
distribution with A = 0.2. 



FIG. 7: Spin excitation spectra in one dimension for (a) 
x = 0.5, (b) x — 0.4, (c) x — 0.3, and (d) x — 0.2, respec- 
tively. The results are for the on-site disorder in the binary 
distribution with A = 0.2. The dashed lines represent the 
Fermi wavenumber at each x. 



features appear more apparently. Especially, in the ID 
case in Fig. [S] (b) , the detailed structure of the anomalies 
can be observed. There, we have anticrossing with gap 
opening and shadow-band like features. From this sys- 
tematic study for different dimensions, we speculate that 
the branching observed in higher dimensions is a remnant 
of this anticrossing. 

The anticrossing with the gap opening reminds us of 
the electron dispersion in systems with spin density wave 
or collective phonon mode. In these systems, magnons 
or phonons scatter the electrons and open a gap with 
the anticrossing in the dispersion. This analogy will be 
discussed further in Sec. IIII D II 



3. Doping dependence 

Here, we study the doping dependence in a system- 
atic way. Figure [5] shows the results for different hole 
concentrations x in 3D cases. We consider here the on- 
site disorder whose distribution is given by the binary 
form of Eq. with A = 0.2. Figure El shows that 
the anomalous features become apparent as x decreases. 
This is probably because the disorder strength becomes 
relatively large compared to the kinetic energy of elec- 
trons which becomes small as x decreases. Moreover, we 
note that the positions of the anomalies appear to shift 
as x. 

In order to observe the x dependence more clearly, we 
study the ID systems under the same conditions. Fig- 
ure shows the results. As indicated in the figures, the 
position of the anomalous features clearly shift as the 
change of x, and is around the Fermi wavenumber k-p 
which is given by tt(1 — x) in this ID case (&f is for the 
spinlcss fermion in the perfectly polarized ground state). 
We note also that the relatively weak anomalies can be 
observed additionally at q ~ n— kp, 2kp, 2(ir— kp), and so 
on. These indicate that the anomalies are closely related 
to the fermionic degrees of freedom. This observation 
provides a key to understand the origin of the anomalies. 
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D. Origin of the anomaly 

1. Friedel oscillation 

We summarize the aspects of the anomalies of the spin 
excitation spectrum which are observed in the previous 
section, (i) There appear additional broadening, branch- 
ing, anticrossing with gap opening, and shadow band, 
(ii) These anomalous features are observed universally 
for different types of the disorder, spatial dimensions and 
doping concentrations, (iii) The anomalies become more 
apparent in lower dimensions, (iv) They appear most 
conspicuously at q ~ kp while the additional anomalies 
can be seen at q ~ Q — kp,2kp,2(Q ~ kp), etc where 
Q = (tt,7v • •)• 

As mentioned in Sec. lIII 0"2l the anticrossing structure 
suggests a scattering mechanism for magnons. Moreover, 
as shown in Sec. IIII C 31 this scattering is characterized 
by the Fermi wavenumber kp, which suggests that the 
scatterers are closely related with the fermionic degrees of 
freedom. From these points, we consider that the origin 
of the anomalies is the Friedel oscillation. 42 

The mechanism is as follows. In the presence of disor- 
der, electrons tend to screen the effects of disorder, which 
leads to the charge density correlation with the wavenum- 
ber 2kp. This is the so-called Friedel oscillation. 42 In this 
DE system, the charge density wave is equivalent to the 
spin density wave since the ground state is the fully polar- 
ized ferromagnetic state (half-metallic stated) . Thus, we 
have the 2kp spin density correlation due to the disorder 
though it is nonmagnetic. The spin density oscillation 
couples to magnon excitations and scatters them. This 
leads to the anticrossing features near q = kp. (Note 
that the scattering wavenumber 2kp corresponds to the 
scattering between q ~ kp and — kp.) The Friedel oscilla- 
tion itself is a general phenomenon in fermionic systems 
with disorder, and its effect becomes relatively weaker 
in higher dimensions. These aspects are consistent with 
our observations mentioned above. Therefore, we con- 
clude that the Friedel oscillation is the origin of the spin 
excitation anomalies in this system. 

2. Single impurity problem 

In order to confirm the picture based on the Friedel 
oscillation, we examine the problem with a single impu- 
rity. In this case, the charge density correlation is clearly 
found around the impurity position. Figure [S| (a) shows 
the density distribution in the real space when we intro- 
duce a single impurity in the ID system with N = 256. 
Here, Uj = (cjcj) and we put the on-site potential at 
i = 128 with relatively large £j = —5 to enlarge the 
impurity effect in the spectrum. The electron density 
oscillates around the impurity site, which results in the 
singularity at q = 2kp in the density correlation function 
as shown in Fig. [S] (b). Note that q — 0.6ir corresponds 
to 2kp with the hole Fermi wavenumber kp — 0.3tt (or 



FIG. 8: (a) Local electron density when the on-site potential 
is finite only at the site i = 128. (b) Density correlation 
function, (c) Spin excitation spectrum. The dashed line in 
(c) represents the Fermi wavenumber. 



2-7T — 2kp with the electron Fermi wavenumber kp = 0.7tt). 
Here, we define the density correlation function as 

In this case of the single impurity also, as shown in Fig. [8] 
(c) , the spin excitation spectrum shows weak but obvious 
anomalies at k = kp,n—kp, etc., which are similar to the 
results in the previous section IIII CI The position of the 
anomalies is solely determined by kp while their strength 
depends on the number and the strength of the disor- 
der. This strongly supports our scenario of the Friedel 
oscillation. 



3. Comparison with random exchange Heisenberg model 

In Sec. Ill Dl we discussed the correspondence between 
the DE model and the Heisenberg spin model within the 
spin wave approximation. There, Eq. I|18|) gives a relation 
between these two models, where the exchange coupling 
Jij is given by the expectation value (ejej). In order to 
show the fermionic effect in this quantity explicitly, we 
study here the spin excitation in the Heisenberg model 
with completely random exchange couplings. Namely, 
we consider the ferromagnetic Heisenberg model with 

(r) 

the nearest- neighbor exchange couplings which take 
random values from bond to bond and no longer satisfy 
Eq. (|18|) . This will elucidate a contrast from the random 
DE model Within the lowest order of 1/5 expansion, 
this model leads to the effective spin-wave Hamiltonian 
in the form 

Hiti = -2Sj\f , H u = 2SJ2 4^ ■ ( 42 ) 

j 

(See Appendix |A*1) Although this effective Hamiltonian 
has a similar structure to the magnon self-energy and 
satisfies the sum rule Eq. I|14[l . it is nothing to do with 
the fermionic degrees of freedom in the DE model. We 
calculate the spin excitation from Eqs. (|42|l in ID. A typ- 

ical spectrum is shown in Fig. |9] Here, we choose J^- 

(r) 

randomly from the range of 1 — A < Jh' < 1 + A with 
A = 0.25. We find no notable anomalies except for the 
broadening for the entire spectrum in contrast to the re- 
sults in Sec. lIII (Jl This result supports the above picture 
and elucidates the relevance of the fermionic degrees of 
freedom of itinerant electrons to the spin excitation in 
the random DE model. 
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FIG. 9: Spin excitation spectrum in the Heisenberg spin 
model with random exchange couplings. See the text for de- 
tails. 
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FIG. 10: Spin stiffness at x — 0.3 as a function of A. Cir- 
cles, squares, and triangles represent the results in one, two, 
and three dimensions, respectively. Open and filled symbols 
correspond to the cases of the bond disorder and the on-site 
disorder, respectively. The lines are guides for the eye. 



E. Quantitative aspects of the spectrum 

As shown in Sec. IIII CI although the spin excitation 
spectrum exhibits significant anomalies near the Fermi 
wavenumber fep, it appears to be always single-peaked 
near the zone center q = 0. This allows us to apply the 
spectral function analysis in Sec. Ill C J 21 in this regime. 
In this section, we discuss the quantitative aspects of 
the excitation spectrum by applying the spectral func- 
tion analysis to the numerical results in comparison with 
the analytical expressions obtained in Sees. Ill C 31 and 
III C 41 Throughout this section, we consider the binary 
distribution i|38|) for both cases of the bond disorder and 
the on-site disorder. 



1. Spin wave stiffness and magnon bandwidth 

As shown in Eq. I|29[l . the stiffness A is proportional 
to the kinetic energy of electrons per site in our analysis. 
Figure plots the stiffness as a function of A at x = 0.3, 
which is calculated from the magnon self-energy by using 
Eq. 126(1 . In all dimensions, the stiffness decreases as the 
disorder becomes strong. The asymptotic dependence in 
the limit of A — > appears to be A cx A 2 . 

Next we consider the bandwidth of the spin excitation. 
In the pure system without disorder, the total bandwidth 
of the spin wave excitation E sw is determined by the spin 
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0.0 0.1 0.2 0.3 

A 

FIG. 11: Total bandwidth of the spin excitation spectrum in 
three dimension with the on-site disorder. Circles, squares, 
and crosses represent the results for x = 0.4, 0.3, and 0.2, 
respectively. The lines are guides for the eye. 

stiffness A in the form E sw = 2zAm& Thus, from Eq. (ZBfr. 
the spin wave bandwidth E sw is proportional to the ki- 
netic energy of electrons per site in the absence of dis- 
order. This relation no longer holds in the presence of 
disorder. We estimate the bandwidth E sw by the peak 
energy of the highest excitation at q = (tt, 7t, it) in 3D 
systems. Figure ITT1 plots the estimates of E sw as a func- 
tion of the strength of the on-site disorder A. Although 
the stiffness always decreases as the disorder strength in- 
creases, the change of the bandwidth E sw depends on the 
hole concentration x: E sw decreases for x 0.3, however, 
E sw increases for x ^ 0.3. At x ~ 0.3, E sw remains al- 
most constant in this range of A. This indicates that in 
our DE system at x ~ 0.3, although the spectrum near 
the zone center becomes narrower due to the decrease of 
the stiffness as the disorder becomes stronger, the total 
bandwidth of the spin excitation is almost unchanged. 
Indeed, this behavior is observed in Fig. [21 



2. Excitation energy and linewidth 

The spectral function analysis in Sec. Ill C~2l concludes 
the excitation energy uj sw cx q 2 and the linewidth 7 cx q 
in the limit of q — > 0. Figure IT2l shows the small q behav- 
iors of lu sw and 7 which are calculated from the numeri- 
cal data by using Eqs. 12211 and l|23|l . All the numerical 
results in different dimensions show consistent q depen- 
dences with the analytical predictions. We also confirm 
that there is the incoherent regime where 7 > uj sw near 
q = 0. These justify the applicability of the spectral 
function analysis which is based on the assumption of 
the single-peaked structure of the spectrum. 

Figure EH plots the values of the coefficient of the q lin- 
ear term, ci in Eq. (|34fl as a function of the stiffness A by 
changing the disorder strength A. As A decreases (A in- 
creases), ci increases almost linearly to the decrease of A. 
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FIG. 12: q dependence of the excitation energy (open circles) 
and the linewidth (filled circles) at x = 0.3 in the presence 
of the on-site disorder in the binary distribution A = 0.2 in 
(a) one, (b) two, and (c) three dimensions, respectively, along 
q = (q, 0,0, ■ • •). The gray solid (dashed) lines are the fits to 

e (a). 
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FIG. 13: q linear coefficient ci in the linewidth as a func- 
tion of the spin stiffness in (a) one, (b) two, and (c) three 
dimensions, respectively. Open and filled symbols correspond 
to the cases of the bond disorder and the on-site disorder, 
respectively. 



The interesting point is that the coefficient ci increases 
much faster in the case of the bond disorder than the 
on-site disorder. This indicates that the bond disorder 
induces larger q linear term, namely, stronger incoher- 
ence in the magnon excitation than the on-site disorder. 

This point is clearly demonstrated when we plot the 
linewidth as a function of the excitation energy as shown 
in Figs. PHI and IT51 For the bond disorder fFig.Hl|. the 

1/2 

behavior 7 oc lo S w dominates the spectrum. On the con- 

___ 1/2 

trary, for the on-site disorder (Fig. H3|) , the w S w part is 
observed only in the small q regime, and the behavior 
7 oc lo sw is dominant in the wide region of q. There, 

we have a crossover from the incoherent behavior with 

1/2 

7 oc w s 4 oc q to the marginally coherent behavior with 



7 oc u> sw oc q 2 . Qualitative features are universal irre- 
spective of the spatial dimensions. 



We find here the qualitatively different behavior be- 
tween the bond disorder and the on-site disorder in the 
crossover from the incoherent to the marginally coherent 
regime. Considering this difference and its origin, we will 
discuss a general tendency of this crossover and a realistic 
type of the disorder in real compounds in Sec. IIV Cl 



11 



y 




0.0 0.1 0.2 0.3 0.4 0.5 0.6 




0.00 0.05 0.10 0.15 0.20 




0.00 0.02 0.04 0.06 0.08 0.10 



FIG. 14: Linewidth as a function of excitation energy in 
the case of the bond disorder. The results are for x = 0.3 
in (a) one, (b) two, and (c) three dimensions, respectively, 
along q — (q, 0, 0, • ■ ■). Circles, squares, and crosses represent 
A = 0.05, 0.1, and 0.15, respectively. 
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FIG. 15: Linewidth as a function of excitation energy in 
the case of the on-site disorder. The results are for x — 0.3 
in (a) one, (b) two, and (c) three dimensions, respectively, 
along q — (q,0,0, ■ ■ ■). Circles, squares, and crosses represent 
A = 0.1,0.2, and 0.3, respectively. 



IV. DISCUSSIONS 
A. Spin excitation anomalies in low-Tc manganites 

In high-Tc manganites, such as Lai-^Sr^MnOa and 
Lai_a;Pba;Mn03 near x — 0.3, the spin excitation spec- 
trum shows a cosine like dispersion^S!^ which can be well 
explained by the pure DE model without disorder (Fig. [21 
(a))i22i Contrast to this canonical DE behavior, low- 
Tc manganites such as Pri-^Sr^MnOs, Lai^^Ca^MnOs, 
and Ndi-^Sr^MnOs show the following anomalous fea- 
tures in the spin excitation spectrum. 

One of the remarkable anomalies is called the soft- 



ening. Hwang et al. reported a significant softening 
of the spin wave dispersion near the zone boundary in 
Pri-^Sr^MnOa at x = 0.37i2£ Later, Dai et al. confirmed 
this result in several materials, and observed the optical 
phonon mode in the energy range where the softening 
occurs^ A similar spectrum with the optical phonon 
branch has been obtained also in the bilayer manganite, 
La2-2xSri + 2 2 ;Mn207 at x = 0.4,-4- j n this compound, 
however, a qualitatively different spectrum, which is a 
cosine like form without any softening, has been observed 
in the same chemical formula but in a different sample 
crystal^ This suggests the strong sample dependence of 
this softening behavior. 
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Another anomaly is the broadeningiSLS&Si In com- 
pounds which show the softening, the linewidth of the 
excitation becomes significantly large in the region where 
the softening is observed. An important aspect is that 
the linewidth remains finite even at the lowest tempera- 
ture where the ferromagnetic moment is fully saturated. 
This point will be discussed in Sec. IIV B 21 

In compounds with relatively small x (low doping), 
the spin excitation shows another anomaly such as the 
gap opening and/or the anticrossing. 28 ' 3 ^ Note that the 
ground state in these compounds is ferromagnetic insu- 
lator and not metal. The phase diagram in the small 
x region shows complicated transitions in structural and 
magnetic degrees of freedom accompanied by a dimen- 
sional crossover 

As shown in the results in Sec. IIIII our simple model 
JU, i.e., the DE model with disorder, reproduces all 
the above anomalies, at least, in a quantitative level. 
The softening appears in the lower energy branch in the 
anticrossing. The additional broadening is observed in 
the large q regime. These anomalies become more con- 
spicuous in lower doped case and in lower dimensions. 
These agreements suggest that the disorder is a candi- 
date to give a reasonable description on the large spec- 
tral changes from the high-Tc to the low- To manganites. 
However, besides our disorder scenario, there have been 
many other theoretical proposals for the present problem. 
We compare our results with those by other scenarios in 
the following. 



B. Comparison with other theoretical results 

1. Softening 

Several mechanisms have been proposed for the soften- 
ing near the zone boundary. One is the purely magnetic 
origin, i.e., due to competitions between the DE ferro- 
magnetic and the superexchange antiferromagnetic inter- 
actions between localized spins including longer ranged 
couplings^ Other scenarios take account of a coupling 
to other degrees of freedom. One proposal is the coupling 
to orbital fluctuations, in which orbital degrees of free- 
dom in the doubly degenerate e g orbitals are explicitly 
taken into account^ Another proposal is the coupling 
to phonons which leads to the anticrossing between the 
phonon excitation and the magnon excitation^! 

Although all these scenarios including our disorder sce- 
nario can reproduce the lower 'softened' energy excita- 
tion in experimental results, there are some qualitative 
differences. In our results, the softening does not oc- 
cur in the strict sense, because the lower energy excita- 
tion appears just as a result of the anticrossing. We still 
have the higher energy branch near the original disper- 
sion in the absence of the disorder. If one follows only 
the lower branch, it looks like the softening as pointed 
out in Sec. IIII C II On the contrary, the mechanisms of 
the magnetic origin and of the magnon-orbital coupling 



origin cause the softening of the original dispersion itself, 
hence there is no higher energy branch. In the scenario 
of the magnon-phonon coupling, the spectrum shows a 
similar structure to our anticrossing. However, in this 
case, the lower excitation near the zone boundary is orig- 
inally the optical phonon mode, which extends over the 
whole Brillouin zone with almost flat dispersion (optical 
mode). There should be a finite energy excitation even 
at the zone center additionally to the gapless magnon ex- 
citation. This is different from our results in which the 
spectrum is single-peaked near the zone center. These 
qualitative differences can be distinguished experimen- 
tally by studying the high energy excitation in the wide 
region of the Brillouin zone. This kind of experimental 
tests will be discussed and summarized in Sec. IIVDI 

Besides the above qualitative differences, our disorder 
scenario is different from the other proposals quantita- 
tively. In other scenarios, it should be noted that the 
relative importance of the additional couplings is deter- 
mined by the change of the DE interaction, namely, the 
bandwidth, since the strength of them does not change 
so much by the ^4-site substitution. In this sense, these 
scenarios attribute the origin of the change of the spec- 
trum from high-Tc to low-Tc compounds to the change 
of the bandwidth. As discussed previously, however, the 
actual change of the bandwidth in these compounds is 
very small. Therefore, it appears to be difficult for these 
scenarios to explain the large changes of the spin exci- 
tations in a quantitative manner. On the contrary, the 
strength of the disorder can be largely affected by the A- 
site substitution even if the change of the averaged lattice 
structure is small. This appears to favor our disorder sce- 
nario. 



2. Broadening 

As discussed in Sec. IIV Al the experimental results 
show the intrinsic linewidth at the lowest temperature. 
This suggests that the spin wave excitation is no longer 
the eigenstate of the system at the ground state. We 
discuss here two different mechanisms for this zero- 
temperature broadening. One is effective even in the 
absence of disorder and the other is our disorder scenario. 

Even in the pure DE model without disorder, the spin 
wave is not the exact eigenstate of the system^ The 
coupling with charge degrees of freedom leads to a fi- 
nite lifetime for the spin wave excitations. This zero- 
temperature broadening in the pure DE model has been 
discussed by the spin wave approximation^^ In the 
lowest order of the 1/S expansion, the DE model can 
be mapped to the pure Heisenberg model as discussed 
in Appendix 1X1 Hence, the spin wave is the exact eigen- 
state and the linewidth is zero in 0(1/ S)~ In the higher 
order of 1/S, however, we have a finite imaginary part in 
the magnon self-energy, which leads to a linewidth in the 
spin excitation spectrum. Up to 0(1/ S 2 ), this linewidth 
is predicted to show the q dependence as 7 oc q d + 3 in 
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the limit of q — > 0, where d is the spatial dimension of 
the system. Thus, the linewidth is predicted to behave 
as 7 oc q e in 3D and 7 oc q 5 in 2D systems. 

In the presence of disorder, our results indicate that 
the spin excitation spectrum shows a finite linewidth 
even in the lowest order of 1/S expansion. As discussed 
111 sec. rnrm this is an inhomogencous broadening and 
distinguished from the above broadening due to the fi- 
nite lifetime of quasiparticle excitations. In our results, 
the linewidth shows the q linear behavior in the limit of 
q — > 0, and crossovers to the q 2 behavior as q increases. 
These behaviors do not depend on the spatial dimension 
d. Therefore, it is possible to determine which scenario is 
relevant in real compounds by studying the q dependence 
of the linewidth in both 3D and 2D systems. 

In 2D (bilayer) compounds, recently, detailed analyses 
have been done on the q dependence of the linewidth. 45 
The marginally coherent behavior 7 oc u> sw oc q 2 is ob- 
served in the wide region of q. It is difficult to explain this 
q dependence by the former mechanism which predicts 
7 oc q b in 2D. In our disorder scenario, the marginally co- 
herent behavior is actually found as shown in Sec. lIII E2l 
This supports our scenario. 

Unfortunately, the linewidth in the small q region has 
not been fully investigated experimentally, especially in 
3D compounds in the ferromagnetic metallic regime^SS 
Detailed information near the zone center especially for 
q <S 7r/5 is lacking. Further experimental study is desired 
to examine the asymptotic q dependence. 

In our results, the crossover between the incoherent 
regime and the marginally coherent regime appears to 
depend on the details of the disorder as mentioned in 
Sec. IIII E 21 The dominance of the marginally coher- 
ent behavior in experimental results provides us a micro- 
scopic picture on the disorder in real compounds. This 
will be discussed in Sec. II V d in detail. 



3. Anticrossing and gap opening 

As mentioned in Sec. II V Al the gap opening and/or 
the anticrossing have been experimentally observed only 
in the low doped region (x ^ 0.2) where Tq is largely re- 
duced from the value at x ~ 0.3i22*22i In this region, the 
kinetic energy becomes small partly because the num- 
ber of carriers decreases. There, the system is insu- 
lating with ferromagnetism or spin canting as well as 
the orthorhombic lattice distortioni^ 1 ^!^! Hence, there 
should be competitions among several interactions such 
as ferromagnetic DE, antiferromagnetic superexchange, 
and electron-lattice interaction. 

We note that, in the 3D system, the values of A in 
the present calculations are much smaller than the crit- 
ical value for the Anderson localizationj 4 ^^ hence, the 
ground state remains the ferromagnetic metal in our re- 
sults. In lower dimensions, of course, electrons are easily 
localized at small A. The ferromagnetic or the spin cant- 
ing state indicates the relevance of the DE interaction 



even in the insulating states. The charge/spin oscillation 
due to the disorder is expected even in these insulators al- 
though the Fermi wavenumber is no longer well defined. 
Hence, we expect similar anomalies in these insulating 
regions. 

Besides our disorder scenario, we have at least two 
mechanisms which cause the anticrossing and/or the gap 
opening in the spin excitation spectrum. One is the su- 
perlattice magnetic structure such as the spin canting. 
In this case, the spin excitation spectrum shows a gap 
at the wavenumber for the corresponding magnetic pe- 
riodicity. The spectrum consists of a single branch in 
the extended Brillouin zone and deviates from the cosine 
like form only near the gap. The small gap observed at 
q = (0, 0, 7r/2) in Lao.ssSro.isMnOa appears to be ex- 
plained by this scenario since the corresponding super- 
lattice peak is observed^ This behavior should be dis- 
tinguished from the anticrossing or the branching in our 
results. 

The other mechanism is the magnon-phonon coupling 
as mentioned in Sec. 11V B II 3 — In this scenario, the anti- 
crossing with the gap opening appears at the wavenum- 
ber where the magnon dispersion crosses the optical 
phonon branch. The excitation spectrum looks similar 
to our results. Possible experimental tests to distinguish 
these scenarios will be discussed in Sec. IIV Dl 

C. Atomic vs. mesoscopic scale disorder 

Our results in Sec . IIII E2l elucidate a difference in the 
spin excitation between the bond disorder and the on- 
site disorder. In the case of the on-site disorder, the 

1/2 

crossover from the incoherent behavior 7 oc w s 4 oc q to 
the marginally coherent one 7 oc tu sw oc q 2 takes place 
in the relatively small q region, and the latter behavior 
is dominant in the wide region of the spectrum. On the 
other hand, in the case of the bond disorder, the inco- 
herent behavior dominates over the whole q region. The 
bond disorder tends to make the spectrum more incoher- 
ent than the on-site disorder in the present model. 

Since the incoherence comes from the local fluctua- 
tions of Jij in Eq. (|TH|) as discussed in Sec. Ill C 51 this 
difference indicates that the bond disorder leads to rela- 
tively larger fluctuations of than the on-site disorder 
in the present model. This is qualitatively understood 
as follows. In Eq. (|18fl . there are two sources of the lo- 
cal fluctuations of J^. One is the transfer integral Uj 

and the other is the expectation value (c\cj). In the case 
of the on-site disorder, tij is taken to be uniform for all 
the nearest neighbor pairs, hence, the local fluctuations 
come only from {c\cj). This expectation value is disor- 
dered in the presence of disorder, however, it has some 
correlations due to the Friedel oscillation as discussed in 
Sec. lIII Dl On the contrary, in the case of the bond disor- 
der, is taken to be completely independent from bond 
to bond in our model. This local fluctuation of £y is di- 
rectly reflected in that of Jij in Eq. (|18|) . Even if there 
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are some spatial correlations in (c]cj), the multiplication 
of tij makes Jij more disordered. Thus, in the present 
model, we expect a larger local fluctuations in the case 
of the bond disorder than the on-site disorder. 

However, this does not mean that the bond disorder 
always makes the spin excitation spectrum more inco- 
herent than the on-site disorder in general. The stronger 
incoherence in the present results is due to the specific 
form of Uj in our model. Even in the case of the bond dis- 
order, the incoherence is possibly suppressed if we have 
some spatial correlations in tij. Instead, the important 
conclusion is that the atomic-scale disorder leads to more 
incoherent excitation than the spatially-correlated or the 
mesoscopic-scale disorder. 

As discussed in Sec. IIV B 21 the experimental results 
in the bilayer compounds show the dominance of the 
marginally coherent regime. From the above conclusion, 
this indicates that in real materials the mesoscopic-scale 
disorder is more relevant than the atomic-scale disorder 
if the disorder plays a major role in the spin excitation 
spectrum. Possible sources of the mesoscopic-scale dis- 
order are a large-scale inhomogeneity or clustering of the 
A ions, twin lattice structure, and grain boundaries. 



D. Experimental tests 

We propose here some experiments to test our disorder 
scenario. 

(i) Sample dependence: Our scenario predicts a corre- 
lation between the sample quality and the anomalies in 
the spin excitation spectrum. Even in the same chemi- 
cal formula, the extent of the anomalies may depend on 
the sample quality. The purity of samples can be mea- 
sured independently, for instance, by the magnitude of 
the residual resistivity. In this test, we note that the A- 
site ordered manganites2i*2SiSi will be useful. In these 
materials, the disorder strength due to the alloying effect 
in the ^4-site ions can be tuned to some extent by the 
careful treatment of the syntheses. It is interesting to 
compare the spin excitation spectrum in the ^4-site or- 
dered compound with that in the ordinary solid-solution 
(^4-site disordered) compound in the same chemical for- 
mula. 

(ii) Comparison with the Fermi surface: The spin exci- 
tation anomalies in our results appear strongly near the 
Fermi wavenumber. It is crucial to compare the posi- 
tion of the anomalies with the information of the Fermi 
surface. The Fermi surface can be independently deter- 
mined, for instance, by the angle-resolved photoemission 
spectroscopy. 

(hi) Doping dependence: Closely related to the above 
test (ii), the doping dependence is also a key experiment. 
In our scenario, by controlling the doping concentration, 
the position of the anomalies shifts due to the change of 
the Fermi surface. Moreover, the anomalies become more 
conspicuous for lower doping concentration. 

(iv) High energy excitation: Our scenario predicts the 



anticrossing instead of the softening. Hence, near the 
zone boundary, there remains the higher energy excita- 
tion above the lower energy branch which looks like to 
be softened. As discussed in Sec. IIVB II this point is 
different from several other theoretical proposals. At the 
same time, it is important to identify the optical phonon 
mode, especially near the zone center, for examining the 
relevance of the magnon-phonon coupling as discussed in 
Sec. IIVB II Therefore, it is crucial to study the higher 
energy region over the whole Brillouin zone. 

(v) q dependence of the linewidth: Our results indicate 
that irrespective of the spatial dimension, the linewidth 
of the spectrum is proportional to q near the zone center 
and shows a crossover to q 2 behavior as q increases. As 
discussed in Sec. IIII E 21 this q dependence appears to 
be consistent with the experimental result in the bilayer 
manganite. In all other scenarios, the linewidth in the 
ground state comes from the interplay between charge 
and spin which is discussed in Sec. IIV B 21 This predicts 
qd+3 behavior near the zone center. It is strongly desired 
to examine the linewidth more quantitatively especially 
in 3D compounds. 



V. SUMMARY AND CONCLUDING REMARKS 

In this paper, we have discussed the disorder effects 
on the spin excitation spectrum in the double exchange 
model. We have applied the spin wave approximation 
in the lowest order of 1/5 expansion. The analytical 
expressions are obtained for the spin stiffness, the exci- 
tation energy, and the linewidth. The most important 
result revealed by the extensive numerical calculations 
is that the spin excitation spectrum shows some anoma- 
lies in the presence of the disorder, such as broadening, 
branching, or anticrossing with gap opening. The ori- 
gin of these anomalies is the Friedel oscillation: Disorder 
causes the 2kp oscillation of the spin and charge density 
in the fully polarized ferromagnetic state, which scatters 
the magnons to cause the anticrossing in their dispersion. 
Thus, the spin excitation anomalies are the consequence 
of the strong interplay between spin and charge degrees 
of freedom. Our results have been compared with ex- 
perimental results in the ^4-site substituted manganitcs 
which shows relatively low-Tc, and the agreement is sat- 
isfactory. We have also compared our results with other 
theoretical proposals for the anomalous spin excitation 
in these manganites, and clarified the advantages of our 
disorder scenario. 

Another important finding is the incoherence of the 
magnon excitation. We have shown that near the zone 
center the linewidth is proportional to the magnitude of 
the wavenumber q, which becomes larger than the exci- 
tation energy which scales to q 2 . This incoherence comes 
from the local fluctuations of the transfer energy of itin- 
erant electrons. As q increases, we have obtained the 
crossover from this incoherent behavior to the marginally 
coherent one in which both the linewidth and the exci- 
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tation energy are proportional to q 2 . We found that this 
crossover behavior depends on the nature of the disor- 
der. Comparing with experimental results in which the 
marginally coherent behavior is dominant, we have con- 
cluded that in real materials the spatially-correlated or 
mesoscopic-scale disorder is relevant compared to the lo- 
cal or atomic-scale disorder. 

All the above anomalous features are obtained in the 
lowest order of the 1/5 expansion, 0(1/5*). Up to 
0(1/ S), we have shown that the double exchange model 
with disorder is effectively mapped to the Heisenberg spin 
model with random exchange couplings. In the absence 
of the disorder, the corresponding Heisenberg model has 
the uniform nearest neighbor exchanges, whose spin ex- 
citation spectrum shows the cosine dispersion. 26,38 Even 
in this pure case, in the higher order of 1/5, there are 
some deviations from this cosine dispersion such as the 
magnon bandwidth narrowing or the softening, and the 
broadeningi 3 ^ 3 ^ 6 " However, these higher order correc- 
tions are known to be small near the hole concentration 
x = 0.3 where we are interested in^i Moreover, the spin 
excitation anomalies we obtained here are the primary ef- 
fects in the 1/S expansion compared to the higher-order 
corrections. Therefore, we believe that our results are 
relevant to understand the anomalous spin excitations in 
real compounds. Our results give a comprehensive under- 
standing of the systematic changes of the spin excitation 
spectrum in the ^4-site substituted manganites. More- 
over, the disorder explains well the experimental facts 
which cannot be understood only by the bandwidth con- 
trol, such as the rapid decrease of Tc~ Therefore, we 
conclude that the disorder plays an important role in the 
j4-site substitution of the ionic radius control in CMR 
manganites. 

We comment on the relevance of another additional 
couplings due to the influence of the multicritical phe- 
nomena. As decreasing the ionic size of the 4-site atoms, 
finally the system becomes the charge-ordered insulator 
with the orbital and lattice orderings, the ferromagnetic 
insulator, the antiferromagnetic insulator, or the spin 
glass insulator Ai 6 ^*^^ The phase diagram often shows 
the multicritical behavior. These phase transitions might 
not be explained by the disorder alone. In the close vicin- 
ity of the multicritical phase boundary, we expect large 
fluctuations and critical enhancement of additional cou- 
plings, which possibly accelerate the suppression of Tq 
and enhance the anomalies in the spin excitation spec- 
trum. However, the multicritical phenomena is essen- 
tially the first order phase transition. Hence, far from 
the phase boundary, for instance, in LCMO compounds, 
we believe that these additional couplings becomes irrel- 
evant and the disorder plays a major role. 
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APPENDIX A 

In this Appendix, we show that the magnon self-energy 
(|llfl is equivalent to the effective spin-wave Hamiltonian 
for the Heisenberg spin model in the lowest order of the 
1/5 expansion^ 

We begin with model Q in the limit of Jn/t — > oo. In 
this limit, since the conduction electron spin <x is com- 
pletely parallel to the localized spin at each site, states 
with a antiparallel to S are projected out. Then, we 
have the effective spinless-fermion model asii 



H=- J2(Uj4c 3 + h.c.) + £i4 c 



(Al) 



where the transfer integral Uj depends on the relative 
angle of localized spins as 



cos £ cos -0- + sin £ sin °J-e~ { ^^ 
2 2 2 2 



. (A2) 



with Sf = 5 sin 9i cos fa , Sf = 5 sin ^ sin #j , and Sf — 
5 cos 9i. Hereafter, we treat the localized spins S as clas- 
sical objects. The transfer integral in Eq. (|A1|) is a com- 
plex variable. The real and imaginary parts are given 
by 



Re U 




(S + S?)(S + S?) 



respectively. Note that Im tij is antisymmetric while 
Re tij is symmetric, as expected from tij = 

Now we apply the spin wave approximation to 
model IjAlfl . We use here the Holstein-Primakoff trans- 
formation 



S - ajaj, 



(.4 + Oi), Sf 



M/^(a|-a<), (A5) 



and take account of the 1/5 expansion up to the leading 
order of 0(1/ S). Substituting Eqs. into Eqs. 
and HA4jl . we obtain 



Rc ti 
tij 
Im ti 



\(a\a 3 + a] ai - of a* - a]aj), (A6) 



1+45 Wo, 
-£g(a\aj - a] ai ). 



(A7) 
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Hence, the effective magnon-electron Hamiltonian up to 
0(1/ S) is given by 



a , — a]a 



a) aj ) 
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x(cjcj +c]a) 

ij 

+ ^4*- (A8) 

i 

Finally, we obtain the effective magnon Hamiltonian by 
tracing out the fermion degrees of freedom in Eq. I|A8() . 
Up to 0(1/ S), the trace can be calculated as the expec- 
tation value for the ground state without any magnon. 
The result is given by 

H ~ ~ ^(4 c j)( a l a j + a]ai ~ a\ai - ajaj), (A9) 

ij 

up to irrelevant constants. Here we use the general re- 
lation (c\cj) — (CjCi) for the perfectly polarized ground 
state without degeneracy. Note that this effective Hamil- 



tonian is derived generally for both the bond disorder 
and the on-site disorder. From the eigenvalues and the 
eigenvectors for Eq. I|A9|) . we obtain the spectral func- 
tion A(q,u>). Hence, comparing with Eqs. © and llUfl . 
the static part of the magnon self-energy H(uj = 0) is 
equivalent to the effective magnon Hamiltonian (|A9(I . 

This effective Hamiltonian <|A9(1 has the same form 
as that for the Heisenberg spin model 7YHois = 
— 2J2i<j JijSi ' Sj within the leading order of the 1/S 
expansion. This correspondence gives the relation be- 
tween the kinetics of electrons in the DE model and the 
exchange coupling Jj» in the Heisenberg model, which 
is given by Eq. l(T%|l . Thus, up to the leading order of 
the 1/S expansion, the static part of the magnon self- 
energy l|ll|) is the same as the effective spin- wave Hamil- 
tonian for the Heisenberg model. In other words, the spin 
excitation spectrum of the DE model is equivalent to 
that of the Heisenberg model with Jy defined by Eq. i|18|) 
in the lowest order of 1/5 expansion. 
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